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I demonstrate that the amplitude of the high-energy scat- 
tering can be factorized in a product of two independent func- 
tional integrals over "fast" and "slow" fields which interact by 
means of Wilson-line operators - gauge factors ordered along 
the straight lines. 
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The starting point of almost every perturbative QCD 
calculation is a factorization formula of some sort. A clas- 
sical example is the factorization of the structure func- 
tions of deep inelastic scattering into coefficient functions 
and parton densities. The form of factorization is dic- 
tated by process kinematics (for a review, see [Q). In 
case of deep inelastic scattering, there are two different 
scales of transverse momentum and it is therefore natural 
to factorize the amplitude in the product of contributions 
of hard and soft parts coming from the regions of small 
and large transverse momenta, respectively. On the con- 
trary, in the case of high-energy (Regge-type) processes, 
all the transverse momenta are of the same order of mag- 
nitude, but colliding particles strongly differ in rapidity. 
Consequently, it is natural to look for factorization in the 
rapidity space. 

The basic result of the paper is that the high-energy 
scattering amplitude can be factorized in a convolution 
of contributions due to "fast" and "slow" fields. To be 
precise, we choose a certain rapidity 770 to be a "rapidity 
divide" and we call fields with 77 > rj Q fast and fields with 
77 < 770 slow where 770 lies in the region between spectator 
rapidity and target rapidity. (The interpretation of this 
fields as fast and slow is literally true only in the rest 
frame of the target but we will use this terminology for 
any frame). 

Our starting point is the operator expansion for high- 
energy scattering [0 where the explicit integration over 
fast fields gives the coefficient functions for the Wilson- 
line operators representing the integrals over slow fields. 
For a 2=>2 particle scattering in Regge limit s m 2 
(where m is a common mass scale for all other momenta 
in the problem t ~ p\ <~ (p'a) 2 ~ P% ~ (p'b) 2 ~ m2 ) we 
have: 



Mpa,Pb Pa, Pb) = ^2 J d 2 xi...d 2 x n (1) 
C i ^(x 1 ,...x n ){ P B\Tr{U il (x 1 )...U in (x n )}\p' B ) 

(As usual, s = (pa + Pb) 2 and t = (j>a — Pa) 2 )- Here 
Xi (i = 1, 2) are the transverse coordinates (orthogonal to 



both £^4 andps) and Ui(x) — (x)^-£^-U(x) where the 
Wilson-line operator U(x) is the gauge link ordered along 
the infinite straight line corresponding to the "rapidity 
divide" 770- Both coefficient functions and matrix ele- 
ments in Eq. ([!]) depend on the 770 but this dependence 
is canceled in the physical amplitude just as the scale 
/i (separating coefficient functions and matrix elements) 
disappears from the final results for structure functions 
in case of usual factorization. Typically, we have the fac- 
tors r~j (<7 In s/m — 770) coming from the "fast" integral 
and the factors ~ g 2 f]o coming from the "slow" integral 
so they combine in a usual log factor g 2 lns/777 2 . In the 
leading log approximation these factors sum up into the 
BFKL pomeron p], (for a review see r ef. ||). Note, 
however, that unlike usual factorization, the expansion 
(|lf) does not have the additional meaning of perturba- 
tive vs nonperturbative separation - both the coefficient 
functions and the matrix elements have perturbative and 
non-perturbative parts. This happens due to the fact 
that the coupling constant in a scattering process is de- 
termined by the scale of transverse momenta. When we 
perform the usual factorization in hard (k± > fj,) and soft 
(k± < fi) momenta, we calculate the coefficient functions 
perturbatively (because a s (k± > p) is small) whereas 
the matrix elements are non-perturbative. Conversely, 
when we factorize the amplitude in rapidity, both fast 
and slow parts have contributions coming from the re- 
gions of large and small k±. In this sense, coefficient 
functions and matrix elements enter the expansion (|l|) 
on equal footing. We could have integrated first over 
slow fields (having the rapidities close to that of pg) and 
the expansion would have the form: 



A{s,t) = Y J J d 2 xi...d 2 x n D h - i "(x u , 
(p A \Tr{U il (x 1 )...U in (x n )}\p' A ) 



(2) 



In this case, the coefficient functions D are the results 
of integration over slow fields and the matrix elements 
of the U operators contain only the large rapidities 77 > 
7/0. The symmetry between Eqs. (1) and (2) calls for 
a factorization formula which would have this symmetry 
between slow and fast fields in explicit form. 

Our goal is to demonstrate that one can combine the 
operator expansions (0) and (||) in the following way: 



A(s,t) 

(pA\U a ^( Xl ) 



A 2 , 



(3) 
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FIG. 1. Structure of the factorization formula. Dashed, 
solid, and wavy lines denote photons, quarks, and gluons, re- 
spectively. Wilson-line operators are denoted by dotted lines 
and the vector n gives the direction of the "rapidity divide" 
between fast and slow fields. 

where Uf = Tr(X a Ui) (\ a are the Gell-Mann matrices). 
It is possible to rewrite this factorization formula in a 
more visual form if we agree that operators U act only 
on states B and B' and introduce the notation Vi for the 
same operator as Ui only acting on the A and A' states: 

A{s,t) = 

( PA \{PB\ex V (ifc?xV ai (x)U?(x)) \p'a)\p'b) (4) 

In a sense, this formula amounts to writing the coeffi- 
cient functions in Eq. (|I]) (or Eq. (0)) as matrix ele- 
ments of Wilson-line operators. (Such an idea was first 
discussed in ref. ||). Eq. (||) (illustrated in Fig. 1) is 
our main result and the rest of the paper is devoted to 
the derivation of this formula and the discussion of its 
possible applications. 

Let us now briefly remind how to obtain the operator 
expansion (Q). For simplicity, consider the classical ex- 
ample of high-energy scattering of virtual photons with 
virtualities ~ — to 2 . 

A(s,t) = -i{0\T{j( PA )j(p' A )j(p B )j(p' B )}\0). (5) 

where j{jp) is the Fourier transform of electromagnetic 
current j^(x) multiplied by some suitable polarization 
e fi (p). At high energies it is convenient to use the Su- 
dakov decomposition: p^ — a p p^ + P p pt^ + p ± where 
Pi and P2 are the light- like vectors close to pa and ps, 
respectively (p A = p^-p^p\/s, p B =p^-p^p B /s). We 
want to integrate over the fields with a > a where a is 
defined in such a way that the corresponding rapidity is 
?7o- (In explicit form ijo = In f where cr = — ). The re- 
sult of the integration will be given by Green functions of 
the fast particles in slow "external" fields ^| (see also ref. 
0). Since the fast particle moves along a straight-line 
classical trajectory, the propagator is proportional to the 
straight-line ordered gauge factor U. For example, when 
x + > 0, y + < it has the form [||: 



G(l , 9) = i ^ i(z ,,_|^^ x) _ti_ <„, 

We use the notations z, = z^p 1 ^ and z* = z^pi^ which 
are essentially identical to the light-front coordinates 
2+ = z^/^/s, z_ = z m /^fs. The Wilson-line operator 
U is defined as 

= [oopi +xx,-oopi +xx] (J) 

where [x,y] is the straight-line ordered gauge link sus- 
pended between the points x and y: 

[x, y]=Pexp (igf* du(x - y^A^ux + (1 - u)y)j (8) 

The origin of Eq. (|^) is more clear in the rest frame of 
the "A" photon. Then the quark is slow and the exter- 
nal fields are approaching this quark at high speed. Due 
to the Lorentz contraction, these fields are squeezed in a 
shock wave located at = 0. Therefore, the propagator 
of the quark in this shock-wave background is a prod- 
uct of three factors which reflect (i) free propagation from 
x to the shock wave (ii) instantaneous interaction with 
the shock wave which is described by the operator U(z±), 
and (iii) free propagation from the point of interaction z 
to the final destination y. 

The propagation of the quark-antiquark pair in the 
shock-wave background is described by the product of 
two propagators of Eq. (^) type which contain two 
Wilson-line factors U(z)U^(z') where z' is the point 
where the antiquark crosses the shock wave. If we sub- 
stitute this quark-antiquark propagator in the original 
expression for the amplitude (|^) we obtain Q : 

/ d 4 zd 4 2e iw - x+i «-*<T{j(o; + z)j{z)}) A 
* f%ihHp±,q±)Tr{U(p±)Ui(q ± -p ± )} (9) 

where U(p±) is the Fourier transform of U(x±) and 
I(p±,q±) is the so-called "impact factor" which is a 
function of p\,p±_ • qx> and photon virtuality ||, [||. 
Thus, we have reproduced the leading term in the ex- 
pansion (To recognize it, note that U(x±)W (y±) — 

Pexp |— ig dz i Ui(zj_)^ where the precise form of the 

path between points x±_ and y± does not matter since 
this is actually a formula for the gauge link in a pure 
gauge field Ui(z±)). 

Note that formally we have obtained the operators 
U ordered along the light-like lines. Matrix elements 
of such operators contain divergent longitudinal integra- 
tions which reflect the fact that light-like gauge factor 
corresponds to a quark moving with speed of light (i.e., 
with infinite energy). As demonstrated in 0, we may 
regularize this divergence by changing the slope of the 
supporting line: if we wish the longitudinal integration 
stop at i] = t]q, we should order our gauge factors U 
along a line parallel to n = ap\ + ap2- Then the co- 
efficient functions in front of Wilson-line operators will 



U(x±) 
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contain logarithms ~ g 2 lnl/a. For example, there are 
corrections of such type to the impact factor I(p, q) and 
if we sum them, the impact factor will be replaced by 

(g 2 In l/cr)" /C"/(p, q) where K, is the BFKL kernel. 

In order to understand how this expansion can be gen- 
erated by the factorization formula of Eq. (||) type we 
have to rederive the operator expansion in axial gauge 
A, = with an additional condition ^4*1^ = ~o = (the 
existence of such a gauge was illustrated in |j] by an ex- 
plicit construction). It is important to note that with 
power accuracy (up to corrections ~ a) our gauge con- 
dition may be replaced by e^A^ = 0. In this gauge the 
coefficient functions are given by Feynman diagrams in 
the external field 



Bi(x) = Ui(x x )Q{x. 



B. = B* = (10) 



which is a gauge rotation of our shock wave (it is easy to 
see that the only nonzero component of the field strength 
tensor F,i(x) — Ui(x±)6(x*) corresponds to shock wave). 
The Green functions in external field ([To]) can be ob- 
tained from a generating functional with a source re- 
sponsible for this external field. Normally, the source for 
given external field A^ is just J v — D^F^ V so in our case 
the only non- vanishing contribution is J*(B) = D l Fi*. 
However, we have a problem because the field which we 
try to create by this source does not decrease at infinity. 
To illustrate the problem, suppose that we use another 
light-like gauge A* = for a calculation of the propa- 
gators in the external field (|lO[). In this case, the only 
would-be nonzero contribution to the source term in the 
functional integral D l Fi t A* vanishes, and it looks like 
we do not need a source at all to generate the field 
(This is of course wrong since B^ is not the classical solu- 
tion). What it really means is that the source in this case 
lies entirely at the infinity. Indeed, when we are trying 
to make an external field A in the functional integral by 
the source J M we need to make a shift A^ — > Ap + A^ in 
the functional integral 

/ VA exp {iS(A) - ifd 4 xJ* {x)A a » (a) } (11) 

after which the linear term D^F fiV A v cancels with our 
source term J^A^ and the terms quadratic in A make 
the Green functions in the external field A. (Note that 
the classical action S(A) for our external field A = B 
( ^p| ) vanishes). However, in order to reduce the linear 
term J d 4 xF^ D^Av in the functional integral to the 
form J d^xD^F^A 1 '(x) we need to make an integration 
by parts, and if the external field does not decrease there 
will be additional surface terms at infinity. In our case 
we are trying to make the external field A = B so the 
linear term which need to be canceled by the source is 

| / dx,dx a ,d 2 x±F im D*A l = / dx*d 2 x±F i ,A l \ x x '~ ^ oo 

It comes entirely from the boundaries of integration. If 
we recall that in our case F,i(x) — Ui(x±_)S(x*) we can 
finally rewrite the linear term as 



/ d 2 x 1 _U i {x 1 _){A i (-oop 2 + x±) - A l (cx>p 2 + x±)} (12) 

The source term which we must add to the exponent in 
the functional integral to cancel the linear term after the 
shift is given by Eq. (|l2|) with the minus sign. Thus, 
Feynman diagrams in the external field ( |l0| ) in the light- 
like gauge A* = are generated by the functional integral 

T>Aexp{iS(A)+ (13) 

tfd 2 x ± U ai (x 1 _)[AU^P2 + x±) - A ai {~^p 2 + ssl)]} 

In an arbitrary gauge the source term in the exponent in 
Eq. ( |l3| ) can be rewritten in the form 

2i J d 2 x±Ti{U l (x±) dv[—oop 2 + x±,vp 2 + x±] 

F„(vp 2 + x ± )[vp 2 + xj_ , -oop 2 + x_l}} (14) 

Thus, we have found the generating functional for our 
Feynman diagrams in the external field (11). However, 
it is easy to see (by inspection of the first rung of BFKL 
ladder diagram) that the longitudinal integrals over a in 
these diagrams will be unrestricted from below while we 
need the restriction a > a. Fortunately, we already faced 
that problem on the other side - in matrix elements of 
operators U and we have solved it by changing the slope 
of the supporting line. Similarly to the case of matrix 
elements, it can be demonstrated that if we want the 
logarithmical integrations over large a to stop at a = a, 
we need to order the gauge factors in Eq.(|l4) along the 
same vector n = ap\ + ap2 , cf. Eq. (||) . Therefore, the 
final form of the generating functional for the Feynman 
diagrams (with a > a cutoff) in the external field (11) is 



VAV^V^e 



iS(A,*)+i J d 2 x ± U"(x ± )V t a (x ± ) 



(15) 



where 



Vi(x±) — f^ ^ dv[—oon + x±, vn + x±] 

n^F^vn + x±)[vn + x_l, —oon + x_l] (16) 

and V® = Tr(A a Vi) as usual. For completeness, we have 
added integration over quark fields so S(A, \&) is the full 
QCD action. 

Now we can assemble the different parts of the factor- 
ization formula (ji|). We have written down the generat- 
ing functional integral for the diagrams with a > a in 
the external fields with a < a and what remains now is 
to write down the integral over these "external" fields. 
Since this integral is completely independent of ( |l5| ) we 
will use a different notation B and \ for the a < a fields. 
We have: 



VADme iS ^j{p A )j{p' A )j(p B )j{p' B ) 



(17) 



/ VAV^e 3 ^ j ( PA )j (p' A ) j VBVxDx 
3{pb)]{p'bY S{B ' x) exp [i f d 2 x ± U a \x ± )V t a (x ± )} 
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The operator K in an arbitrary gauge is given by the 
same formula (|l6| ) as operator Vi with the only difference 
that the gauge links and F,i are constructed from the 
fields This is our main result (jj) in the functional 
integral representation. 

The functional integrals over A fields give logarithms 
of the type g 2 In 1 /a while the integrals over slow B fields 
give powers of g 2 \n(as/m 2 ). With logarithmic accuracy, 
they add up to g 2 In s/m 2 . However, there will be addi- 
tional terms ~ g 2 due to mismatch coming from the re- 
gion of integration near the dividing point a ~ a where 
the details of the cutoff in the matrix elements of the 
operators U and V become important. Therefore, one 
should expect the corrections of order of g 2 to the effec- 
tive action J dx±U l Vi. 

In conclusion let us discuss possible uses of this ap- 
proach. First thing which comes to mind is to use the 
factorization formula for the analysis of high-energy ef- 
fective action. Consider another rapidity rj' in the re- 
gion between 770 and In m 2 / s. If we use the factorization 
formula ( |l7| ) once more, this time dividing between the 
rapidities greater and smaller than r]' , we get: 



VAe^j(p A )j(p' A )j(p B )j(p' B 



labels 



In 



(x—x')' 



In 
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(x-x') 2 {y-y')' 2 



where we we use the notation V^{x) = g§r^°(a;) etc. 
The first term here looks like the corresponding term in 
the factorization formula (|l7|) - only the directions of the 
supporting lines are now strongly different. The second 
term is the first-order expression for the reggeization of 
the gluon and the third term is the two-reggeon Lipa- 
tov's Hamiltonian [[12j responsible for BFKL logarithms. 

This approach can be also used for the study of high- 
energy heavy-ion collisions since it is well suited for the 
study of the interaction of two colliding shock waves. In- 
deed, for heavy-ion collisions the coupling constant may 
be relatively small due to high density (see |il|] ) . On the 
other hand, the fields produced by colliding ions are large 
so that the product gA is not small - which means that 
the Wilson- line gauge factors V and Y are of order of I . 
In this case we need to know not only a couple of the 
first few terms in the expansion of the effective action, 
but the whole series. 

(18) Acknowledgements. The author is grateful to L.N. 
Lipatov and A.V. Radyushkin for valuable discussions. 
\VAe} siyA) 7 (va )i(v' ) VBe iS( - B "' j(nn)j(n' )e iSait ^ Vi ' Yi '^ This work was supported by the US Department of En- 
J ' J ergy under contract DE-AC05-84ER40150. 

where the effective action for the rapidity interval be- 
tween rj and rf is defined as 



(19) 



J J d 2 x ± V ai (x ± )Ui(x±)+i J d 2 x ± W ai (x ± )Y i a (x±) 



(For brevity, we do not display the quark fields). In this 
formula the operators [/, are constructed from C fields 
while the operators Wi (made from C fields) and Yi (made 
from B fields) are given by the same Eq. (|n|) with gauge 
links aligned along the direction n' = u'pi + &'p2 corre- 
sponding to the rapidity 77' (as usual, lncr'/cr' = 77' where 
a' = m 2 /so'). 

The formula (|l^) gives a rigorous definition of the ef- 
fective action for a given interval in rapidity (cf. ref. 
|^|). Next step would be to perform the integrations over 
the longitudinal momenta in the r.h.s. of Eq. jl9| ) and 
obtain the answer for the integration over our rapidity 
region (from 77 to 77') in terms of two-dimensional the- 
ory in the transverse coordinate space which hopefully 
would give us the unitarization of the BFKL pomeron. 
At present, it is not known how to do this. For illustra- 
tion, let us present a couple of first terms in the effective 
action El, El: 
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